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We give an analytic construction of a class of two-qubit gate pulse sequences that act on five of the 
six spin-1 particles used to encode a pair of exchange-only three-spin qubits. Within this class, the 
problem of gate construction reduces to that of finding a smaller sequence that acts on four spins and 
is subject to a simple constraint. The optimal sequence satisfying this constraint yields a two-qubit 
gate sequence equivalent to that found numerically by Fong and Wandzura. Our construction is 
sufficiently simple that it can be carried out entirely with pen, paper, and knowledge of a few basic 
facts about quantum spin. We thereby analytically derive the Fong-Wandzura sequence that has so 
far escaped intuitive explanation. 


Control of the Heisenberg exchange coupling between 
pairs of spin-i particles is a useful resource for carrying 
out quantum gates in a quantum computer [T] . If qubits 
are encoded in the Hilbert space of one [T] or two [5] 
spin-i particles, resources beyond exchange are required 
for universal quantum computation. However, controlled 
exchange alone is universal if each qubit is encoded in 
the Hilbert space of at least three spin-i particles [3H5]. 

Semiconductor quantum dots with trapped electrons 
are promising systems for manipulating spin-f particles 
[5]. Controlled exchange between pairs of electron spins 
in quantum dots has been demonstrated [7], and used 
to manipulate a variety of three-spin encoded qubits i- 
[13], including the so-called resonant exchange qubit ni- 
dS] which maintains qubit encoding by keeping the ex¬ 
change interaction “always on” within each qubit (see 
also El)- Here we focus on the case of exchange-only 
quantum computation where the exchange interaction is 
kept completely off except when being pulsed, i.e. adi- 
abatically switched on and off, between pairs of spins. 
It is then necessary to design pulse sequences that carry 
out quantum gates on encoded qubits without resulting 
in leakage out of the encoded qubit space OdUdiHia. 

To assess any quantum computation scheme one ulti¬ 
mately needs to know the minimal cost of carrying out 
quantum gates. For exchange-only quantum computa¬ 
tion using pulse sequences, single-qubit gate sequences 
are theoretically understood but there is little true un¬ 
derstanding regarding optimization of two-qubit gate se¬ 
quences. The main difficulty comes from the constraint 
of no leakage combined with the large search space of uni¬ 
tary operators acting on the six spins encoding a pair of 
three-spin qubits. Not surprisingly, the shortest known 
pulse sequence for an entangling two-qubit gate due to 
Fong and Wandzura [20] has been found by a numerical 
search algorithm which offers little insight into its deriva¬ 
tion. Furthermore, existing analytic derivations of less 
optimal sequences are lengthy and complicated 

In this Letter we present an analytic construction of a 
class of pulse sequences that carry out two-qubit gates for 
exchange-only quantum computation. These sequences 
are built out of smaller sequences that act on only four 



FIG. 1. (a) Three-spin encoding for logical qubit states with 
total spin I, |a) with a = 0,1, and noncomputational state 
with total spin |, |nc). (b) Two qubits in states a and b and a 
diagram highlighting the five rightmost spins with total spin 
/=|or|. 


spins and satisfy a certain constraint. We show that when 
the most efficient of these smaller sequences is used the 
result is equivalent to the Fong-Wandzura sequence. Our 
guiding principle throughout is to avoid as much as possi¬ 
ble complicated calculations and use only the most basic 
facts about quantum spin [23] . 

Because we only consider the action of rotationally in¬ 
variant operators, it is sufficient to describe quantum 
states of multiple spins using only total spin quantum 
numbers. Accordingly, we employ the notation used in 
[22] in which each spin-^ particle is represented by the 
symbol • and groups of particles are enclosed in ovals la¬ 
beled by the total spin of the enclosed particles. Figure 
l^a) shows the three-spin qubit encoding of [5] in this 
notation. In the text we write these states using paren¬ 
theses instead of ovals, so the computational qubit states 
shown in Fig. [^a) are = |a) where a = 0 and 

I define the standard basis. Figure [^a) also shows the 
non-computational state (• • •) 3/2 = \nc). Figure 0b) 
shows six spins encoding two qubits and highlights the 
five spins the pulse sequences we consider act on, where 
the total spin / can be either ^ or |. 

Consider the exchange Hamiltonian H = JS^ ■ Sj act¬ 
ing on two spin-f particles (••)a whose Hilbert space is 
spanned by the two states with total spin a = 0 and 
I. (Because only total spin quantum numbers are rele¬ 
vant in our analysis, we treat, for example, the three-fold 
degenerate a = I state as a single state.) Pulsing this 
Hamiltonian for a duration t measured in units of l/(7rJ) 
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FIG. 2. (a) Elementary exchange pulse of duration r = 0 or 1, 
referred to as an r-pulse in the text, and matrix representation 
of the resulting operation in the basis a = {0,1}, where m = 1 
or —1 for r = 0 or 1, respectively, (b) Sequence of three r- 
pulses and two explicit SWAP pulses (r-pulses with r = 1) 
acting on three spins and the matrix representation of the 
resulting operation in the basis ac = {0|,l||l|}. 

results in the time evolution operator (up to an irrelevant 
overall phase factor), 


FIG. 3. (Color online) (a) Operation, R, acting on a spin-| 
particle, •, and an effective spin-| particle, ★, which repre¬ 
sents the three spins encoding the target qubit. The matrix 
representation of this operation, which generalizes the r-pulse 
of Fig.[^a), is given in the effective basis d = {0,1} where 1 is 
the 2x2 identity, and M is a 2 x 2 matrix where = 1, both 
of which act on the target qubit contained in ★. (b) Pulse 
sequence which generalizes that shown in Fig. Sb). The ma¬ 
trix representation of the resulting operation is given in the 
effective basis af = {0|, 


17y(f) = diag(I,e “‘), (I) 

where the matrix representation is given in the a = {0,1} 
basis. We take t G [0, 2) for which the inverse pulse has 
duration 2 — t. 

Two exchange pulses square to the identity and play 
a fundamental role in our construction. The durations 
of these pulses, which we denote r, can only be either 0 
or 1, and we refer to them as r-pulses. For r = 0 an r- 
pulse is simply the identity, while for r = 1 it is a SWAP 
operation, which is equivalent to physically exchanging 
two spins [24]. Figure [^a) shows an r-pulse acting on 
two spins in the state (••)a using standard notation with 
exchange pulses represented by double arrows labeled by 
duration. The corresponding matrix representation of 
the resulting operation, also given in Fig. [^a), shows 
that applying an r-pulse multiplies the a = 0 state by 
1 and the a = 1 state by m, where m = lorm = —1 
for r = 0 or r = 1, respectively. In both cases = 1, 
reflecting the fact that the r-pulses square to the identity. 

A pulse sequence that acts on three spins and consists 
of three r-pulses (with either r = 0 or 1) acting on the 
bottom two spins and two explicit SWAP pulses acting 
on the top two spins is shown in Fig. [^b). For r = 0 
the two explicit SWAP pulses square to the identity. For 
r = 1 the sequence consists of five SWAP pulses which, 
when viewed as spin permutations, are readily seen to 
be equivalent to a single SWAP pulse acting on the top 
two spins. In both cases, the effect of the sequence is 
to multiply the state by 1 if a = 0, and m if 

a = 1, regardless of the value of c, where, as in Fig. [^a), 
TO = 1 or TO = —1 for r = 0 or r = 1, respectively. 
The corresponding matrix representation is also given in 
Fig.§b). 

We seek pulse sequences which act on the five spins 
((••)a(*(**)6)i/2)/ highlighted in Fig. j^b) and carry 
out leakage-free two-qubit gates. For reasons that will 


become clear we refer to the qubit with state label a 
as the control qubit and the qubit with state label b as 
the target qubit. Our construction is based on using a 
smaller sequence which acts on the four rightmost spins 
in Fig. 0b)- This smaller sequence carries out an opera¬ 
tion we denote R which, as seen shortly, is closely related 
to an r-pulse. One requirement we place on R is that it 
not result in any leakage of the target qubit into its non- 
computational state. We can therefore work within an 
effective Hilbert space in which the three spins encoding 
the target qubit are replaced by a single effective spin-| 
particle, 

(•(••)ft)i/2 —t ★. (2) 

Matrix elements of operations acting on any collection 
of spins including ^ are then elevated from numbers to 
2x2 blocks that act on the Hilbert space of the target 
qubit hidden within 

We require that when R is applied to the state (•★)d 
it act on the target qubit with the identity 1 if d = 0, 
and a matrix M, with = 1, if d = 1, as also shown 
in the corresponding matrix representation of R given in 
Fig.j^a). Such R operations can be viewed as generalized 
r-pulses where the matrix elements 1 and to, with to^ = 
1, of Fig. [^a) have been elevated to the 2x2 matrices 1 
and M, with = 1 in Fig. 0a). 

This view of R as an elevated r-pulse suggests the five- 
pulse sequence of Fig. 0b) can also be elevated to the 
sequence shown in Fig.^b). This sequence acts on the 
effective Hilbert space spanned by the states ((••)a-A-)/ 
with a/ = 0^, and l| and consists of three R opera¬ 
tions acting on the central spin and the effective spin ★ 
and two SWAP pulses acting on the top two spins. The 
only 2x2 block element in the matrix representation of R 
which is not proportional to the identity is M. Because 
= 1, when evaluating the matrix representation for 
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the full sequence, each block matrix element must be of 
the form ag 1 + ai M. The coefficients ag and ai for 
each block element are completely determined by the two 
cases M = ±1, which are equivalent to the cases m = ±1 
in Fig.[|;b). It follows that the matrix representation of 
the operation carried out by this sequence in the effec¬ 
tive af = {0|, l^llH basis is that given in Fig. 3 b), i.e. 
an elevated version of the matrix shown in Fig. 2 'b). To 
prove this result we have only used the fact that = 1. 
It therefore holds not just for the straightforward cases 
M = ±1, but also for any matrix of the form M = h ■ cr 
where n is a real-valued unit vector and cr = ( ctj ,, ay,az) 
is the Pauli vector. 

The pulse sequence shown in Fig. [^b) acting on the 
two qubits of Fig.[^b) applies the identity 1 to the target 
qubit when the state of the control qubit is a = 0, and 
applies the matrix M to the target qubit when the state 
of the control qubit is a = 1, regardless of the value of 
/. The matrix representation of the operation carried 
out by this sequence can then be given in the standard 
two-qubit basis ab = {00,01,10,11} as 


U 2 qubit = diag{l,M). (3) 

For M = ±1 the resulting gate is not entangling. How¬ 
ever, for M = n • (T the sequence enacts a leakage-free 
controlled-(n -cr) gate which is equivalent to a controlled- 
NOT (CNOT) gate (for which h = x), up to single-qubit 
rotations. 

Abandoning the notation we now consider R act¬ 
ing on the four-spin Hilbert space spanned by the states 
where, since c is initially d can only be 
either 0 or 1. The requirements on R needed to con¬ 
struct a controlled-(n-cr) gate are that it must i) preserve 
the quantum number c, and ii) in the restricted Hilbert 
space with c = ^, have the form shown in Fig. j^a) with 
M = h - cr. 

To construct a sequence for R we introduce a new op¬ 
eration V which satisfies the constraint 

(((••)i(**)i)i|^^l(*(* • •)3/2)i) = 0 (4) 

depicted in Fig. Qa). As shown below, inserting any 
V satisfying into the sequence shown in Fig. ib) 
results in an R operation with M = n • cr. Letting 
Uij{t) denote an exchange pulse of duration t acting on 
spins i and j, as defined in Q, the sequence for R can 
be written F“^{ 7 i 2 (l)t/ 34 ( 1 )F, using the spin labeling 
of Fig. Sb). The matrix representation of the central 
two SWAP pulses 1 / 12 ( 1 ) 1 / 34 ( 1 ) in the ((••)b'(••);,)(; ba¬ 
sis with state ordering bb'd = { 000 , 110 | 011 , 101 , 111 } is, 

t/i 2 (l)// 34 (l) = diag(l, 1| - 1, -1,1). (5) 

In the d = 0 sector 1 / 12 ( 1 ) 1 / 34 ( 1 ) acts as the identity, and 
thus R also acts as the identity since V and V~^ cancel 
one another. In the d = 1 sector, Q implies that V maps 
the c = I state (•(•••) 3 / 2 )i entirely into the b'b = 01,10 



FIG. 4. (Color online) Constructing a pulse sequence for R. 
(a) Constraint which must be satisfied by an operation, V, 
used in the construction, (b) Sequence in which V, its inverse, 
and two SWAP pulses are used to carry out an R operation, 
(c) Evaluation of the matrix element ([^ for the case tit 2 = 11, 
as described in the text. (All quantum numbers are the same 
as those in (a), but are omitted for readability.) (d) One of 
two two-pulse solutions of (a) for V which can be used to 
construct R using (b). 


subspace. The two central SWAP pulses then apply a 
phase factor of —1 to any state in this subspace, and so 
after applying V~^ the net effect of the full sequence will 
be to multiply the c = | state by —1. The fact that R 
maps the c = | state onto itself immediately implies that 
R also maps the c = / subspace onto itself, and thus leads 
to no leakage of the target qubit. Furthermore, since the 
trace of 1 / 12 ( 1 ) 1 / 34 ( 1 ) in the d = 1 sector is —1 (see (H)), 
the trace of the full sequence 1 / 1 / 12 ( 1 ) 1 / 34 ( 1 ) 1 / ^ in this 
sector is also —1. Thus, since the c = | matrix element of 
the full sequence is —1, the trace of the operation acting 
on the c = / subspace must be 0. Finally, because the 
two SWAP pulses 1 / 12 ( 1 ) 1 / 34 ( 1 ) square to the identity the 
full sequence for R squares to the identity. The operation 
carried out by this sequence on the c = | subspace in the 
d = 1 sector must therefore also square to the identity 
and, because it is traceless, it must have the form M = 
n • cr. 

The set of pulse sequences V that satisfy the constraint 
Q can be used to construct an infinite class of sequences 
resulting in two-qubit gates which are locally equivalent 
to CNOT. We now show that the fewest number of pulses 
needed to satisfy Q is two and for this optimal case 
the resulting two-qubit gate sequence is equivalent to the 
Fong-Wandzura sequence. 

Without loss of generality we take V = 1 / 23 ( 12 ) 1 / 12 ( 11 ) 
|25j . While the values of li and I 2 for which V satisfies 
can be found by brute force calculation, we determine 
them here by a simple “pen and paper” procedure. This 
procedure is based on the observation that Q implies 

(((••)l(**)l)l|^23(l2)l/l2(ll)|(»(» • •)3/2)l) = 

a -b + 7e-“‘= -b (6) 

where the coefficients a, /3, 7 , and S can be found by 
evaluating the four simple cases lil 2 = 00, 01,10, and 11. 
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FIG. 5. (Color online) Sequence of Fig.j^b) acting on two en¬ 
coded qubits resulting in a leakage-free controlled-(n-cr) gate. 
Also shown is the full sequence obtained by first inserting the 
optimal sequence for V from Fig.|^d) into the sequence for R 
from Fig. Sb) and inserting the result into the sequence from 
Fig.§b). This sequence is equivalent to the Fong-Wandzura 
sequence. 

For the case tit 2 = 00 the matrix element (§ is simply 
equal to the overlap F = (((••)i(**)i)i|(*(* • •) 3 / 2 )i) 
[2S] . For both cases tit 2 = 01 and 10 there is a single 
SWAP pulse which can be applied either directly to the 
left (for tit 2 = 01 ) or right (for tit 2 = 10 ) four-spin 
state in Fig. |^a). Since each SWAP pulse then acts on 
a pair of spins with total spin 1 , the result is an overall 
— 1 using the phase convention of ([^. In both cases the 
matrix element § is thus equal to — F. 

The only non-trivial case is tit 2 = 11 for which both 
pulses are SWAP pulses. A method for evaluating (§ in 
this case is sketched in Fig. |^c). First, a pair of SWAP 
pulses which combine to the identity, 1 / 24 ( 1 ) 1724 ( 1 ) = 1 , 
is inserted at the start of the sequence. We then view 
the four SWAP pulses as physical particle exchanges. It 
is irrelevant that the effect of particle exchange differs 
from that of a SWAP pulse by a factor of —1 because 
there are an even number of SWAP pulses. Applying 
one of the exchanges acting on spins 2 and 4 to the state 
(•(• • •) 3 / 2 )i then gives a factor -|- 1 , since the two spins 
being exchanged have total spin 1. The remaining three 
exchanges can then be applied to the state ((••)i(**)i)i 
where, referring to the Fig. j^c), they result in a permu¬ 
tation which exchanges the bottom two spins (red oval) 
with the top two spins (black oval). The net effect is 
therefore to exchange two spin -1 objects with total spin 
1 and, as a result, the right state acquires a factor of — 1 . 
Thus the tit 2 = 11 matrix element ([^ is equal to — F. 

Having evaluated the left-hand side of (|^ for the four 
cases tit 2 = 00 , 01 , 10 , and 11 , the coefficients appearing 
on the right-hand side are easily found to be —a = j3 = 
7 = d = F/2. For these coefficients there are only two 
solutions which satisfy Q, tit 2 = 5 §,§ 5 - Figure |^d) 
shows the resulting sequence for the first solution, which 
consists of one -s/SWAP (t = 1) and one inverse VSWAP 
(t = |) pulse. 

Figure shows the pulse sequence obtained by insert¬ 
ing the sequence for V from Fig. id) into Fig. ®b) 
and inserting the resulting sequence for R into Fig. [^b). 
This sequence can be applied to a linear array of spins 
with nearest-neighbor pulses and carries out a controlled- 
(n-cr) gate consisting of eight SWAP pulses, six VSWAP 


pulses, and six inverse \/SWAP pulses [37]. After single¬ 
qubit rotations are removed, the core Fong-Wandzura se¬ 
quence as published in [2^ consists of six SWAP pulses, 
three -s/SWAP pulses, and nine inverse Vs WAP pulses. 
Our sequence can be converted into the Fong-Wandzura 
sequence through a series of elementary manipulations in 
which pairs of SWAP pulses are inserted (as in Fig. 0c)) 
or removed from the sequence, and single SWAP pulses 
are pulled past other pulses and in some cases combined 
with VSWAP pulses to form inverse Vs WAP pulses (and 
vice versa). The same manipulations can be used to pro¬ 
duce sequences applicable to spin geometries other than 
linear [21]. These manipulations preserve the fact that 
there are twelve nontrivial (i.e. not SWAP) pulses in 
these sequences as well as the parity of the sum of the 
number of SWAP and VSWAP pulses. This parity is 
odd for the Fong-Wandzura sequence and even for our 
construction, consistent with the fact that a single-qubit 
operation corresponding to a single SWAP pulse must be 
added to the core Fong-Wandzura sequence to produce a 
controlled-(fi • cr) gate. 

In summary, we have analytically constructed a class 
of pulse sequences for carrying out two-qubit gates for 
exchange-only quantum computation. These sequences 
can be viewed as elevated versions of the simple three- 
spin sequences shown in Fig. [^b) which consist entirely 
of SWAP operations. To carry out this elevation we in¬ 
troduced the four-spin sequence R which is itself built out 
of a smaller sequence V which satisfies the constraint Q . 
When the shortest pulse sequence for V is plugged back 
into the full two-qubit sequence the result is equivalent 
to the Fong-Wandzura sequence. 
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